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Abstract — Since the performance of the teleoperation systems 
can be considerably degraded by time-delay of communication 
channels and uncertainty in various parts of such systems, 
the main objectives of the controller design in loads of different 
structures of the bilateral teleoperation system are to preserve 
stability and tracking performance of these systems in spite 
of aforementioned sources of uncertainty. In this paper, a 
new robust PID controller will be designed based on H„ control 
theory by using the Linear Matrix Inequality (LMI) approach. 
Therefore, the problem of a Robust PID controller design can 
be regarded as a special case of the output-feedback controller 
via employing some sorts of changes in control and system 
parameters. To show the effectiveness of the proposed 
controller, the robust PID controller is compared with the 
multiobjective H/H„ one. The main feature of the suggested 
structure is its ability to control the teleoperation system via 
using the simplest structure in which two signals will be 
transmitted to control the teleoperation system. In addition, 
use of PID controller has more practical applications in 
industrial units, due to its simplicity in implementation and 
capability to predict the time responses caused by changes in 
control parameters. 

Index Terms — Teleoperation systems, Robust LMI based PR) 
controller, Multiobjective H/H„, Optimization. 

I. Introduction 

There has been a considerable growth in use of 
teleoperation systems in different areas such as: telesurgery, 
telemanipulation, space mission, nuclear power station, under 
sea research, Etc. A typical bilateral teleoperation system 
depicted in fig. 1 consists of five important parts: a human 
operator, a master manipulator, communication channels, a 
slave manipulator and remote environment. The master 
manipulator is directly drawn by the human operator and its 
position/velocity is transmitted to the slave site via 
communication channel. The slave manipulator has to track 
the position/velocity of master one and communicate 
reactions of the remote (task) environment to the master 
manipulator as the reflected force. 

Although using the bilateral structure in teleoperation 
systems increases ability of the human operator to control 
the teleoperation system, if there was a long distance between 
slave and master sites, performance and stability of the 
overall system would be deeply affected. The performance 
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of a teleoperation system is addressed by an index called 
transparency. Transparency is defined as a match between 
forces that is exerted by the human operator and one which is 
reflected from the task environment. In other words, 
transparency is a capability of a teleoperation system to 
represent unchanged dynamics of the remote environment to 
the human operator [ 1 ] . Because of the existing uncertainties 
in dynamics of the system and time-delay in communication 
channels a compromise is required to be ensured between 
transparency and stability [2]. Since an increase in stability 
margins of the system leads to a decrease in transparency, 
control of a bilateral teleoperation system requires a delicate 
trade-off between two foregoing requirements namely 
transparency and robust stability. 

Until now several control schemes have been applied to 
bilateral teleoperation systems so next some of these schemes 
are considered briefly. H,, control theory was used by Sano et 
al. in 2000 [3]. They used four sensors in order to measure 
positions and forces of master and slave systems. In 2002, a 
new control method consisting of smith predictor and wave 
variables was proposed by Ganjefar et al. [4] for optimizing 
the performance of the teleportation system against the large 
variable time-delay. An optimal H 2 procedure was proposed 
for teleportation systems in 2004 by Boukhnifer and Ferreira 
[5] . Sunny et al. [6] designed an adaptive position and force 
stabilizing controller for bilateral teleportation system in 2005. 
In the same year Ganjefar and Miri [7] used optimization 
control method in teleportation systems. One year after that, 
ShaSadeghi et al. [8] proposed a new control structure based 
on adaptive inverse method while at this year Sirouspoor 
and Shahdi [9] employed predictor controller for teleportation 
systems. Tavakoli et al. suggested a modeling and stability 
analysis for discrete teleoperation system in 2008 [10]. While 
a sliding mode bilateral control was proposed by Moreau et 
al. for Master-Slave pneumatic servo systems [11]. Novel 
adaptive-based methods have been dedicated to various 
structures of teleoperation systems [12-15] to reduce the 
destructive effects of time delay and uncertainty on 
performance of these systems. Additionally, different 
disturbance based methods have been proposed for bilateral 
teleoperation systems to preserve the stability and 
performance [16-18]. 
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Figure 1. General structure of the bilateral teleoperation system 

Recently, LMI-based approaches have been employed to 
deal with stability and performance problems [19]. In 2008, an 
LMI approach was proposed to design robust H„ and L con- 
trollers for bilateral teleoperation systems by M. ShaSadeghi 
et al. [20]. A teleoperation system with asymmetric time-de- 
lays was also studied in the form of LMI [21] in 2009. The 
stability condition based on LMI has been used to optimize 
the maximum permitted value of time-delay. The major short- 
age of the above papers is that the complete transparency 
has not been achieved. 

This paper presents a novel method to design a robust 
PID controller based on an LMI approach to control a bilat- 
eral teleoperation system with uncertainty in large time-de- 
lay and parameters of the slave and remote environment. The 
main objectives of the controller design are to achieve com- 
plete transparency and robust stability simultaneously for 
the closed-loop system. The proposed method consists of 
two local controllers. One controller is placed on the remote 
site called slave controller and has to be responsible for track- 
ing the commands of the master system. Referring to the 
uncertainty in remote site, a new robust PID controller is 
supposed to be used in the remote site. The second LMI 
based robust PID controller which is called master controller 
is positioned on the local site to be in charge of transparency 
and stability of the Closed-loop. In order to show the effec- 
tiveness of the designed controller the simulation results will 
be compared with the conventional multiobjective HJH 
controller. 

The paper is organized as follows: modeling of the 
teleoperation system and time-delay of communication chan 
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nel will be introduced in section 2. Section 3 is devoted to the 
design of the local controllers. Moreover, the stability analy- 
sis of the proposed structure is described at this section. 

The si mul ati on resul tswi 1 1 be shown i n secti on 4. And f i nal I y, 
section 5 draws conclusions and some suggestions for the 
futurework must begiven. 

1 1 . Modeling Of The Teleoperation System 

Fig. 2 shows an especial structure [22] of the bilateral 
teleoperation system that uses direct force-measurement at 
the remote site. In this Figure G, C and y denote the transfer 
function of local systems, local controllers, and positions of 
the local manipulators, respectively. Moreover, the subscripts 
m and s is designated to components of the master and slave 
sites, respectively. T and T denote the forward and the 
backward time-delay, whereas K p , K are position and force 
scaling factors. Z is the impedance of the task environment. 
In addition, F ' F , F , F are the human operator, task 
environment and input forces applied to the slave and master 
by their local controllers, and also F is the reflected force 
from task environment. Finally, v represents the sensor noise 
of the force measurement in the remote site. It is worth noting 
that due to use of direct force measurement in considered 
control structure, force sensors are demanded to be used. 
Hence, the noise of measurement sensor in the remote site 
should be considered. This force feedback from remote 
environment makes it easier to control a typical bilateral 
teleoperation system through providing a better perception 
and insight about the remote environment for the human 
operator. This feedback channel, however, may simply lead 
to instability of teleoperation system with regard to increase 
in amount of time-delay. 

A. Modeling of the manipulators 

One degree of freedom teleoperation system has been 
used for modeling of bilateral teleoperation systems in several 
literatures [23]. So dynamics of the master and slave 
manipulators can be described by the state space 
representation of (1). 
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Figure 2. Special Structure of teleoperation system 

vcACEEE 



Long Paper 



Int. J. on Recent Trends in Engineering and Technology, Vol. 9, No. 1 , July 20 1 3 



In (1), M, B are inertia and damping coefficients of the 
manipulators and j=s, m is designated to the slave and master 
respectively. In state space representation of the robots x } 
and x 3 are the velocities of the master and slave manipulators 
whereas x 2 , x 4 are the positions of the master and slave 
manipulators. Additionally,^ and/| are the input forces of the 
controllers that will be exerted to the manipulators in order to 
make desired changes in positions/velocities of them. 
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B. Modeling of the communication time-delay 

It is worth regarding that the forward T and backward 
T time-delays are assumed to be bounded and not necessary 
identical. By assuming T= T +T , 7 can be represented by 
the signal-flow model [24], [25] as follows: 
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Where and 9=(T -T )/(2T ) are real parameters such 
that -ld"MT and Od'W'0.5. Moreover, T and T are the 

' max mm 

upper and lower bound for the time-delay respectively. When 
P varies between -1 to 1, the time-delay takes the values 
between its minimum and maximum values and if T =0, 0=0.5. 

mm 

By using the Laplace transformation, the time-delay in (2) 
can be shown in frequency domain and approximated by the 
first-order Pade of (3). 
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In order to represent an uncertain system into the 
multiplicative form, suppose that the real system belongs to 
a family of plants 77 and can be defined by using the following 
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perturbation to the nominal plant G g : 

G(s) = (l+A(s)W m (s))G (s), VG(*) e n (4) 

In this equation, WJs) is a stable transfer function 
introducing the upper bound of the uncertainty and "(s) 
indicates the admissible uncertainty block, with "(s) „<1. 
Achieving the W (s) involves some numerical searches in 
the frequency domain so that (5) is supported. 



G{ja>) 
Go {Jo}) 
G(jco) 



G (jo) 



l = A(jco)W m (ja>) 
\<\W m (ja>)\ , 



(5) 



Hence, in order to utilize the robust control method, the 
second uncertain part in (3) can be expressed by the 
multiplicative uncertainty structure as (6). 
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Therefore, it is equal to (7.2). 
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This formulation of the communication channel time-delay 
is depicted in fig. 3 and used in H,, framework. 

III. Controller Design 

As it is mentioned before, the main goal of the suggested 
control scheme is to achieve transparency and stability of 
the system versus foregoing sources of instability. This is 
done by designing two local controllers; one is placed in the 
remote site (slave controller) denoted by C and the other 
one is sat in the local site (master controller) represented by 
C . The remote controller guarantees the position/velocity 
tracking while the local controller must preserve force tracking 
(transparency) and stability of the overall system. 

In order to control the bilateral teleoperation system, new 
robust LMI-based PID controllers are established based on 
H„ theory in both sites. In the suggested method, the PID 
controller is designed as a special case of the output-feedback 
controller and then convex optimization will be employed to 
find optimal parameters of the controllers for selected 
weighting functions. A plant with its related weighting 
function can be shown as the generalized structure depicted 
in fig. 4. The state space of the depicted plant is defined by 
(8). 

In (8), u is the control input, w is a vector of exogenous 
signals (reference input signal, disturbance and sensor noise 
signals) and e is the measurement output signals whilst z is a 
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vector of output signals z 2 
performance of the total system. 

x = Ax+ B w w+ Bu 
z = C z x + D zw w + D,u 
e - Cx + D u ,w 



z) T ) related to the 



(8) 



The dynamical output-feedback controller will be 
illustrated with a state space realization given in (9) such that 
makes desired changes in the output vector. 



(9) 
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u = C K C, + D k e 

Furthermore, the state space realization of the closed 
loop can be introduced by (10). 
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From bounded real lemma [26] A2 is stable and H,, norm 
of the desired channel between w and z. in the augmented 
plant is smaller than y if and only if there was a symmetric P 
such that ( 1 1 ) is met. 
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Where B h ., C. and D are related to the desired input w. 
and output z, from the augmented plant in fig. 4. By employing 
the congruence transformation in the form of diag{TI p I,I) [27] 
in which Ilj must be introduced by (12). The none-convex 
problem in (1 1) can be converted to the convex problem with 
new parameters introduced in (13) via exerting congruence 
transformation of ( 1 2) and will be applied in theorem 1 . 
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In (12), A' and M are nonsingular matrices that must 
satisfy M^I-XY, 
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Figure 3. Schematic of the uncertain time-delay 



A k := 
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D k :=D k 

In the next subsections, design of the local controllers 
will be described based on the foregoing principles. 

A. LMI based robust PID controller design 

The considered local systems with their controllers have 
to be stated as the structure in fig. 5. 

In fig. 5 W . is the uncertainty weighting function due to 
variations in the dynamics of the slave manipulator and remote 
environment or change in time-delay which can be achieved 
from (5). W p ., W u . are the performance and controller weighting 
functions and will be chosen based on the objectives of the 
controller design in both sites. In order to design the 
performance weighting function the characteristics of the 
desired output in both frequency or time domain (like cut off 
frequency or settling time or maximum overshoot) should be 
constructed in the frequency domain and so the desired 
transfer function of the system will be established from these 
characteristics [28]. Concerning the nominal performance 
condition W i> .{j(o)Sijco)y. <1 from robust control theory, W . 
will be achieved. Also in order to design W . by considering 
the bounded effect on the controller signal namely 
W {jco)f.(jco)x <1, W . can be constructed from the same method 

as W .. Where i=m, s referring to the master and slave sites. 

pi o 

Then for i=s: 



G s '(s)=G s (s)/{l + Z e G s (s)) 



(14) 



G(s) and Z are introduced before. In fact, equation (14) 
introduces the transfer function of Y/F in fig. 2, and must be 
employed in design of C that is meant to support the position 
tracking of the slave system. Moreover, for i=m: 



G m'0)= G mO) 



C s (s)G s (s)Z e 



l + G s {s)Z e + C s (s)G s {s) 



T ni S ) (15) 



This comes from the Closed-loop formation of the slave 
site with its designed controller. As the master controller is 
considered to guarantee force tracking as well as stability of 
the overall system, the transfer function of F/F must be 
used to meet these objectives. Therefore, by using the 
structure of fig. 3 instead of communication time-delay, (15) 
has to be used in addition to W _ which is connected with 
variation of time-delay. 

So the procedure of controller design for this structure 
includes two regular step: at the first step a controller is 
designed for the slave site (i=s) with related weighting 
functions and then by transmitting the slave Closed-loop to 
the local site the master controller will be established (i=m). 

The structure depicted in fig. 5 could be shown as the 
augmented plant illustrated in fig. 4. The conventional LMI 
based H„ output-feedback controller will be achieved from 
theorem 1 [27]. 

Theoreml: Position tracking for the slave and master 
systems against the uncertainty in parameters of the slave 
system, remote environment, and time-delay could be 
supported if there were X and Y for the following minimization 
LMI. 
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The LMIs of (16) are achieved from (11) by using 
congruence transformation introduced in (12). 

In (16), 77 =77/, and A, B ,C,D are introduced before, 
whereas A k , B k , C k represent transformed parameters of the 
controller of (13) which are used in order to preserve 
linearization of the problem. The real parameters of the 
controller can be constructed from the relations in (17). 
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So next, in order to design a robust PID controller, the 
PID controller design problem should be put into the output- 
feedback controller structure and then the extended theorem2 
can be used to establish the robust PID controller. A PID 
controller signal is founded by ( 1 8). 



u(t) = (K p K, K D )[e \e e'J 



(18) 



Converting e into e ={e +"e e2 ) T 'm the augmented plant, 
the new PID-augmented plant shown in fig. 6 is established 
and the problem of PID controller design can be shown as a 
special case of the output-feedback controller. So the state 
space of the augmented plant is transformed into (19). 



x = Ax + B w w + Bu 
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Where xB is achieved from (20). 
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So by forming the new augmented plant from (21) and 
(23) the transfer function between an especial input w. and 
output z. can be represented by (24). 
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Although A,, B : , C, will be achieved for the new structure 

k k k 

of the PID-augmented plant from therom2, by considering 
the applied structure and regarding that D = (k p k { k D ), the 
coefficient of D k will be enough to guarantee stability and 
performance of the system. 

Regarding the new parameters, an extended version of 
theorem 1 can be used to design a robust PID controller. 

Theorem 2. A robust PID controller could support the 
tracking performance (either position or force) and stability 
of the system if there were ^ , y for the following convex 
LMIs. 

mininmize y 
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In (25), A. -A J' and other parameters are introduced in (26). 
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Proof. The new congruence transformation can be 
introduced by (27) 
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NP and MP are nonsingular matrices with proper 
dimension that must satisfy MPNP T =I-XPYP 

Now applying diag{0 p I,I) to (1 1) will end up in (28). 
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Now considering (24) and (27) theoem2 will be achieved 
so that the relation between the real parameters of the 
controllers and the converted ones are represented by (29). 
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C:=C k M' 
D:= D, 
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Guaranteeing the system stability requires an additional 
constraint D,>0 which is added to the constraints of the 

k 

theoreml. As D k is equal to the coefficient of the PID 
controller this constraint supports the positivity of the 
parameters of the controller and therefore, stabilization of 
the designed controller 




Figure 5. Slave and master sites with weighting functions 

B. LMI based multiobjective HJH,, controller design 

Standard structure of fig. 5 can be shown as the 
generalized plant depicted in fig. 4, and from the generalized 
plant a conventional multiobjective H./H,. controller can be 
constructed from theorem3. 
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Figure 6. Schematic of the PID-augmented plant 
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Theorem3. By considering uncertainty in the 
communication channels and accompanied noise with 
reflected force (received from the remote environment via the 
backward communication channel), force tracking 
(transparency) can be achieved by using a multiobjective 
H/H.. controller in the master site. So a multiobjective H/H.. 
controller with tracking error bound y exists if there were X 
and Y that would satisfy the minimization linear matrix 
inequality in theoreml with following additional constraints: 



<0, 



Table I. Parameters Of System 
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Tr{Q)<v, D +D D k D =0 



Where 77„, 77 2; , IJ 3] , IJ 22 , IJ 32 are introduced before and: 



M n :=C z X+D z C k 



M 32 := C z 



(31) 



And the real parameters of the Controller can be constructed 
from (17). 

Proof. By considering the introduced state space 
representation of the system, H 2 norm of the system will be 
smaller than v if there was a P which satisfies LMIs in (32) 



r A' T P + PA' T PB.'^ 



D„ =0, 



-/ 



-<0, 



c' Q 

trace (Q)<v 



>0, 



(32) 



In (32), Q is an auxiliary parameter applied to maintain 
linearization. By using the transformation of (12) it can be 
converted to (30) and will be added to the constraints of 
theoreml to design a multiobjective /////. controller. 

IV. Simulation Results 

In order to show the effectiveness of the proposed method 
the simulation result of a teleoperation system with parameters 
of Table. 1 will be analyzed. 

Considering the changes of the slave parameters in table. 1 
the multiplicative uncertainty weighting function for the slave 
manipulator and time-delay of communication channels can 
be obtained from (5). 

W ms (s) = 18.05(j+1.035)/(s+5.499) 

W mm (j) = W mT (s) = 1.55/(0.43295 + 1) 
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(33) 
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BJN/m) 


11 
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3 




0.1 


e 


0.5 


m; 
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z t 


1 
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3 




3 


b; 


27 


K p , k, 


1,1 



The bode diagrams of fig .7 show that the uncertainty 
weighting functions resulted from the changes in time-delay 
of communication channels and parameters of the slave site, 
cover the modified range of variations in parameters. 

So with notion the uncertainty weighting function of the 
slave site for the generalized plant with W u =0.0008 (to restrict 
the maximum controller output to below 3), and W =5/ 
(s+0.001 ) (in order to have the slave position without steady 
state error and fast time response), the robust PID controller 
(34) has been achieved from extended theorem2 by using 
numerical computations of the mincx solver in Matlab toolbox 
that support the tracking performance. 



C s (t) = (33.9381 5.8247 5.898xl0" 4 )£ s (?) 
E s (t) = (e s (t) \e s {r)dr e' s (t)J 



(34) 



By employing shown structure in fig. 5 and W =0.2 (in 
order to restrict controller input bellow 5), W m =l/(s+10 ~ 5 ) 
(to impose the desired time response to the output of the 
master manipulator with lower settling time and without 
steady state error), the generalized structure of fig. 6 has 
been used in theoreml in order to reach the tracking 
performance with robustness against time-delay in 
communication channel. 



C m (t) = (0.9277 2.5221xl0~ 5 0.793l)£„,(?) 
E a {t) = {e m {t) \e m (T)dr ej (tjf 



(35) 



Where ep)=yJt-TJ-y s (t) and e m (t)=f h (t)-f e (t-TJ. On the 
other hand by using theorem2 and applying Normalized 
Coprime Factorization (NCF) method, the first 8th order 
multiobjective H/H.. controller will be reduced to the following 
4 th order controller for the master site introduced by (36). 



-3.923 xl(T 5 


-0.2186 


0.0065 


0.0119 








-1.64 


3.625 











-3.625 


-1.64 


1.907 














-5.903 


8 


0.0238 


-3.645 


1.094 


1.996 






(36) 
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Time response of the teleoperation system with designed 
controllers for nominal values of the parameters and time- 
delay are shown at fig. 8 to fig. 12. 

Simulation results show the responses of the system for 
the input signals that is usually used in such systems. As it 
can be seen from fig. 8 to fig. 12 the designed robust PID 
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controller has a faster time response rather than the 
multiobjective H/H,. controller while the control effort is 
acceptable. On the other hand, the PID controller cannot 
reject measurement noise that imposes more control effort to 
this controller. 

Since the simulation results show the responses of the 
system for the nominal value of the parameters u analysis 
has been applied to test the stability and performance of the 
system for different values of the parameters. 

The last two figures of this section indicate the stability 
and performance statuses based on the u analysis. Since the 
values of the u bound related to the stability of the system is 
smaller than one and for performance is near one, the designed 
controllers preserve the global stability regardless the 
uncertainty in dynamics of the slave manipulator and variable 
time-delay in communication channel while the performance 
of the system do not change seriously. 

BaJfpkl tr/thr uncurtain nibm arrj weighting fiincfia 
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Figure 7. The bode diagram of the W and G/G -l and bode diagram 
of the W and approximation of time-delay 
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Figure 8b. Measurement signal (e (f)) 
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Figure 9a. Position of master manipulators {y (()) 




Figure 9b. Position of slave manipulators (y (t)) 
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Figure 9c. Position error signals (e (t)) 
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Figure 8a. Human operator and Reflected forces 
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Figure 10a. Master controller signal (f (£)). 
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Figure 10b. Slave controller signal (J (/)) 
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Figure 11c. Position of slave manipulators (y (t)) 
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10 ■ : 

Figure 12a. i bound of Closed-loop stability and performance with 
designed controller 

Conclusions 

Due to the transmission time-delay in communication 
channel, it is fairly difficult to ensure the Closed-loop stability 
of a teleoperation system, so major concerns of the controller 
design for a bilateral teleoperation system are preserving the 
stability and transparency of the teleoperation system despite 
the time-delay in communication channels. In this paper a 
new robust PID controller has been proposed by applying a 
Linear Matrix Inequality approach in order to guarantee 
complete transparency against the variable time-delay with 
known upper bound and the uncertainty in parameters of the 
slave manipulator and remote environment. Two separated 
controllers have been established to ensure the complete 
transparency and stability of the Closed-loop. In order to 
show the effectiveness of the proposed method, the 
simulation results for both multiobjective H^H„ and robust 
H„ PID controllers in master site have been shown, while the 
slave controller was a robust PID controller in both situations. 
It can be seen from the simulation results that the robust PID 
controller shows acceptable results in comparison with the 
multiobjective H/H,, controller. Furthermore, the PID 
controller has more practical applications. Although the 
designed PID controller cannot reduce the effects of the 
measurement noise, in many practical cases the manipulators 
have the low frequency characteristics, consequently it does 
not effect on the final response of the system intensely. In 
the following researches the noise reduction can be added to 
the robust PID controller by using some stochastic 
calculations to reduce the control efforts of such controllers. 
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